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Abstract 

We estimate the number of solutions of certain diagonal congru- 
ences involving factorials. We use these results to bound exponential 
sums with products of two factorials n\m\ and also derive asymp- 
totic formulas for the number of solutions of various congruences with 
factorials. For example, we prove that the products of two factori- 
als n\m\ with max{n, m} < p l / 2+£ are uniformly distributed modulo 
p, and that any residue class modulo p is representable in the form 
mini + n\ ! + . . . + 77.49! with max{m, n, n\, . . . , 7749} < p 8775 / 8794+e . 

2000 Mathematics Subject Classification: 11A07, 11B65, 11L40. 
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1 Introduction 



Throughout this paper, p is an odd prime. Very little seems to be known, 
or even conjectured, about the distribution of n\ modulo p. In Fll in jH], 
it is conjectured that about p/e of the residue classes a (mod p) are missed 
by the sequence n\. If this were so, the sequence n\ modulo p should assume 
about (1 — l/e)p distinct values. Some results of this spirit have appeared 
in 0. 

The scarcity of heuristic results is probably due to the hardness of com- 
puting factorials. The best know algorithm to compute n! over 2 or modulo p 
takes about n 1 / 2 arithmetic operations in the corresponding ring, see |TJI2|- It 
has been shown in that the complexity of computing factorials is related 
to such deep conjectures of the complexity theory as the algebraic version of 
the P = NP question, see also a nice discussion in [2|. 

Sums of multiplicative characters and various additive and multiplicative 
congruences with factorials have been considered in P HH] • 

In particular, it is has been shown in 7j that for any nonprincipal char- 
acter x modulo p we have 

L+N 

£ X (n\) = O (iV^y^logp) 1 / 4 ) (1) 

n=L+l 

and also that the number of solutions Ie(L, N) of the congruence 

nxl . . . n e \ = n e+ i\ . . . n 2 e\ (modp), L < m, . . . , n 2 e < L + N, (2) 
satisfies the bound 

I e (L,N) <A^-i+2-< (3) 

(provided 0<L<L + N<p). 

Using (Q) and (jSJ), it has been shown in [7j that for any fixed e > 
the products of three factorials ni!n 2 !n 3 !, with max{ni, 
are uniformly distributed modulo p. Here, we obtain an upper bound for 
the additive analogue of (J2J) and use it to estimate double exponential sums 
with products of two factorials. Namely, for integers a, K, L, M and iV we 
consider double exponential sums 

K+M L+N 

W a (K,M;L,N) = $^e(am!n!), 

m=K+l n=L+l 
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where we define 

e(z) = exp(2niz/p). 

In turn, our bound of exponential sums W a (K, M; L, N) lead us to a substan- 
tial improvement of the aforementioned result, showing that the products of 
two factorials n\m\, with max{n, m} = 0(p 1//2+e ), are uniformly distributed 
modulo p. We then combine our new bounds and the bounds (JIJ and (jSJ) to 
study various congruences involving factorials. 
Studying single exponential sums 

L+N 

S a (L,N)= Yl e ( an! ) 

n=L+l 

is of great interest too. Although we have not been able to obtain "individ- 
ual" bounds for these sums, we obtain various bounds "on average", which 
also play a crucial role in our arguments. 

Some results and techniques of [7j have found their applications to study- 
ing prime divisors of n\ + f(n) for various functions /, see |171 ITSj . In 
particular, in they have led to an improvement of a result of Erdos and 
Stewart [S]. We expect that the results of this work will also find some 
applications to various arithmetic questions. 

Throughout the paper, the implied constants in symbols t O > and '<C' may 
occasionally, where obvious, depend on integer parameters k, £, r and a small 
real parameter e > 0, and are absolute otherwise (we recall that f/<V and 
A = 0(B) are both equivalent to the inequality \U\ < cV with some constant 
c> 0). 

Acknowledgements. During the preparation of this paper, F. L. was 
supported in part by grants SEP-CONACYT 37259-E and 37260-E, and I. S. 
was supported in part by ARC grant DP0211459. 
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2 Bounds on the Number of Solutions of Ad- 
ditive Congruences with Factorials 

For integers t > 1, A, L and N with 0<L<L + iV<pwe denote by 
Je(L, N; A) the number of solutions to the congruence 

2^ rij! = raj + A (mod p), L + 1 < n\, . . . , n^t < L + N. 

i=i i=e+i 

We also put J e {L, N) = J t (L, N, 0). 

Our treatment of Je(L, N; A) is based on exponential sums. Accordingly, 
we recall the identity 

y>M = ( ' ! mo ^' (4) 

^-^ I p, if u = U (mod p), 

a=0 v 

which we will repeatedly use, in particular to relate the number of solutions 
of various congruences and exponential sums. 

Theorem 1. Let L and N be integers with 0<L<L + N<p. Then for 
any positive integer £ > 1, the inequality 

j,(AiV;A)«iV 2 ^ 1+1/(m) 

holds. 

Proof. The identity (@J) implies that 

J e (L, N;\) = -J2 \3 a {L, N)\ 2e e (-oA) . (5) 

y a=0 

In particular 

1 

Je(L, N; A) < J e (L,N) = - V \S a (L, N)\ 2i . 
For any integer k > we have 

L+AT 

S a (L,N)= e(a(n + k)\) + 0(k). 

n=L+l 
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Therefore, for any integer K > 0, 

1 K-\ L+N 

S a (L,N) = -J2 E e(a(n + fc)!) + 0(K) 



fc=0 n=L+l 



2/ 



n=L+l k=l \ i=l 

= # £ E e U»!ll(n + i)) +0(4 

n=L+l fc=0 V i=l 

Using the Holder inequality, we derive 

p-1 A"-l L+Af p-1 

j2\s a (L, N r«K-^N^ e E E e (^ ! ^^w)+^ : 

a=0 k\,...,k2t=0 n=L+l a=0 

where 

*jw«(*) = EII("+o- E II( n+< )- 

u=l i=l u=i+l i=l 

The sum over a vanishes, unless 

n\^k u ...,k M (n) = (modp), (6) 

in which case it equals p. 

It is easy to see that ®ki,...,k 2e (X) is a nonconstant polynomial of degree 
O(K), unless (ki, . . . , ke) is a permutation of (ke + i, . . . , kit)-, which happens 
for 0(K e ) choices of < ki, . . . , k^i < K— 1. If $^ ^^(X) is a nonconstant 
polynomial, then © is satisfied for at most K values of n, otherwise we use 
the trivial bound N on the number of solutions in n. 

Because n\ ^ (mod p) for 0<L<n<L + N<p, the total number 
of solutions of ^mL+l<n<L+N and < k\ , . . . , kit < K — 1 is 
0(NK e + K 2£+1 ). 

Thus 

p-i 

J2\Sa(L,N)\ 2£ « K- 2e N 2e - l (NK e + K 2e+l )p + K 2£ p 

a=0 

= (N^K-e + N^K + K^p. 

Taking K = |^jV 1 '^ +1 )J and remarking that with this value of K the last 
term never dominates, we finish the proof. □ 
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Corollary 2. Let 8i = ±1 for each i = 1, . . . , k. Then for any integer X, the 
number of solutions of the congruence 

k 

S y^ j 5,^. = X (mod p), L + 1 < ni, . . . , n k < L + N, 

i=l 

is O (iv fc - 1 + 1 /2ft+i)+i/2(fc 2 +i)) j where kl = \_ k / 2 \, k 2 = |_(Jfe + l)/2j. 

Proof Note that k — k\ + k 2 . By the identity the number J of solutions 
of the above congruence can be expressed via exponential sums as 

. L+N p-1 / / k N 

J = ~ H ^eL ^W-A 

ni,...,n fc =L+l a=0 \ \i=l 



= -^e(-aA)n^ ia (AiV). 

P a=0 i=l 

Since IS^^L, N)\ = \S a (L,N)\, we see that 
p— i p— i 



J <-J2\ S ^ N )\ k = -£|S a (£,iV)| fel |S a (L,iV)| fc2 • 



P S=0 P 6=0 



Using the Cauchy inequality and Theorem [T] (see also ©), we finish the 
proof. □ 



Let Fi(K, M] L, N) denote the number of solutions of the congruence 

l 21 

^^n.;!mj! = njlmj (mod p), 

i=l i=£+l 

K + 1 < mi, . . . , m 2 £ < X + M, L + 1 < m, . . . , n 2 i < L + N. 



(7) 



The condition N 2 > M > N 1 ^ 2 , requested in our next result can be 
substantially relaxed. However, because we are mainly interested in the 
"diagonal case" M = N (for which this condition is always satisfied) and in 
order to avoid some technical complications, we use this condition. 

Theorem 3. Let K, L, M and N be integers with 0<K<K + M<p and 
0<L<L + N<p. For any positive integer £ > 1, such that N 2 > M > 
N 1 / 2 , the following bound holds 

F e (K,M;L,N) « m 2 *' 1 ^ 2 * N 2 ^ 2 ^ . 
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Proof. First of all we note that if we prove the above inequality for N 2 > M > 
N then we are done. Indeed, if N > M > N 1 ' 2 then we have M 2 > N > M, 
and the statement follows from the inequality 

So, let N 2 > M > N. We set H = [M 1-1 / 2 * jVVa(<+i)J . We see that 
N i-i/2t(t+i) ^ H ^ Juf i-i/2«(/+i)_ Theil; by the identity ©, 



M; L,N) 



1 

^E 



2i 



K+M L+N 

e(om!n!) 

m=lf+l n=L+l 
H L+JV 

e(om!n! 

r=l A"+(r-l)A//.ff<m<A _ +rA:/'/.H' n=L+l 



2/ 



Applying the Holder inequality, we obtain 

F £ (K,M;L,N) 



H p-1 



r=l P a=0 



L+V 



e(am!n! 

K"+(r-l)M/iJ<m<.fir+rM/.ff n=L+l 



2C 



Therefore, 



F e (K, M; L, AT) < tf^Q 



where Q is the number of solutions of the congruence J7J) with the additional 
condition that |mj — m 3 -| < M/H, 1 < z < j ' < 2£. Without loss of generality, 
we may suppose that mi = minjmj | 1 < i < 2£}. We denote m\ — m and 
put 



rrii = m + Si 



2<i<2£. 



Then the congruence (J2J) in the new variables takes the form 



2(: 



Hi 



(8) 



i=2 



i=£+l 



where /(to, t) = (m + 1) . . . (m + t) for an integer t > 1 and /(n, 0) = 1. 
The number of solutions of the congruence (jHJ) is collected from two sets 



of variables m and rii, Si, 1 < i < 2£: 
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(i) the first set is such that the left hand side of (JSj) is a polynomial of m 
of degree greater than zero (but less than M/H); 

(ii) the second set consists of those for which the left hand side of (jSJ) is 
constant clS db polynomial of m. 

The number of solutions Q\ of (JBJ) corresponding to the first set is at most 



Qi < N M [j^ + 1 J Jj < (MN/H) . 

For the second set of variables, we have that as a polynomial, the left 
hand side of ((HI) is a constant. Let us numerate s 2 , . . . , s 2 e in an increasing 
order. Then, instead of equation (JHJ) we consider the equation 



21 

n\\ + 5i^2f(m,ri)ni\ = (modp), 

i=2 



with 5j = ±1, and such that 

= n < . . . < r 2 t < M/H. 

Moreover, for each positive integer k < 11 and positive integers e\, . . . , e^, 
we consider solutions with 



(9) 



— V\ — . . . — r ei < T ei _|_i — ... — r ei _|_ e2 < . . . 

< r ei +...+e k _ 1+ i = ...= r ei+ ... +£k < M/H. 

In this case, the vanishing of the polynomial in m on the left hand side of (0) 
leads to the conditions 

nil + 5 2 n 2 \ . . . + 5 ei n ei \ = (modp), 
5 ei+ in ei+1 \ + . . . + 5 ei+e2 n ei+e2 \ = (modp), 

S ei+ ... +ek _ 1+1 n ei+ ... +ek _ 1+1 \ + ... + 5 ei+ ... +efc n ei+ ... +e J = (mod p). 

Certainly, for each solution to the system of congruences (jlOj) there are at 
most M possible values for m. We also note that from (jlOj) it follows that 
ei > 2 for 1 < i < k. In particular, k < £. 
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For each k- dimensional vector e = (ei, . . . , e&) of positive integers such 
that e\ + . . . + efc = 2£, there are O ((M/H) k ~ 1 ) possible integer vectors 
(ri, . . . , r 2 e) satisfying 0. For each such fixed vector (7*1, ... , r 2 e), the number 
of solutions of the system of congruences (fTUj) , by Corollary is at most 

Q ^JJ i Yr,-l+l/2(Le I /2j+l)+l/2(L(e l +l)/2j+l)j = q ^rfU-k+K^ ? 

where 

k 

«(e) = £ 



. , v 2(Le,/2j+l) 2(L(e i + l)/2j+l) 
Therefore, 

Q 2 < max(M/if)^ 1 MAr 2£ - fc+K(e) , 

e 

where the maximum is taken over all integers 1 < k < I and fc-dimensional 
vectors e = (ej, . . . , e&) of integers > 2, 1 < i < k, with e\ + . . . + = 2£. 
If £ = 1 then fc = 1, «(e) = 1/2 and 

Q 2 < MiV 3/2 . 

Therefore, in this case we have 

/•', ( /\ . M : L, N) <<:; M 2 N 2 H ; + MN 3/2 H, 

and the required estimate follows from the choice of H. 
Now, we suppose that i > 2. If = 1, then 



i+l 

If > 2, then trivially 



Hence, 



e < -. 

; ~ 2 



Q 2 « MN^-W+V + iV 2 ^ max {M/HN l ' 2 ) k . 

2<k<2V 



One verifies that for our choice of H and under the condition N e+1 > 
M (whih is always satisfied for I > 2 and iV 2 > M), we have M/HN 1/2 < 1, 
so the term corresponding to k = 2 dominates. Therefore, 

Q 2 « MN 2e - 1 ^ + M 2 N 2i - x H-\ 
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Thus, putting everything together we obtain, 

F e (K,M;L,N) < H 21 ' 1 (Q 1 + Q 2 ) 

« H 2 "- 1 {[MN/Hf 1 + MN 2£ ~ e W + V + M 2 N 2 "- l H- 1 

Since I > 2 then M < jV 3£ /( m ) for N 2 > M. Therefore the first term always 
dominates the third one, and we derive 

F t (K, M; L, N) < M 2i N 2l H~ l + MN 2e ~ e ^ H n ~ l . 

Recalling our choice of H, we finish the proof. □ 

3 Bounds of Double Exponential Sums with 
Factorials 

Unfortunately we are not able to estimate single sums S a (L,N), however 
we obtain nontrivial bounds for double exponential sums with factorials. 
We follow some ideas of Karatsuba [HI HI] and Korobov [HI El|, see also 
Lemma 4 in and the follow-up discussion. 



Theorem 4. Let K, L, M and N be integers with Q<K<K+M<p and 
0<L<L + N<p. Then for any integers k, i > 1, the inequality 



max 

gcd(a,p)=l 



\W a (K,M;L,N)\ « j y 1 - 1 /2^(fc+i)jyi-i/2fc(£+i)pi/2M 



holds. 

Proof. Let Ge(L, N; A) denote the number of solutions to the congruence 

l 

2_. n ^ = A (mod p), L + 1 < ni, . . . , < L + N. 



i=i 



Clearly, 



p— i p— i 

^2G t (L,N;\) = N e and ^G e (L, N; X) 2 = J e (L, N). (11) 

A=0 A=0 



10 



By the Holder inequality, we have 

K+M 

\W a (K,M;L,N)\ e < M^ 1 

m=K+l 



L+N 

E 

n=L+l 



e{am\n\ , 



K+M 



L+N 



= M l ~ l Yl e(am\(n 1 \ + ...+n e \)) 

m=K+l ni,...,ri2i=L+l 
K+M L+N 

= M e ~ 1 d ™ Yl e(am\(n 1 \ + ... + n £ \)) 

m=K+l ni,...,ri2i=L+l 

for some complex numbers d m with |?9 m | — 1, K + l<m<K + M. 
Therefore, 

K+M p-l 

\W a {K,M-L,N)\^M 1 - 1 Y,G t (L,N',\)0 m e(a\m\). 

m=K+l A=0 

Applying the Holder inequality again, we derive 

p-l K+M 



\W a (K,M; L,N)\ 2ki < M 2 ^ 1 ^ 



A=0 



m=K+l 



2k 



'p-l 



2k-l 



2fe/(2fe-l) 



vA=0 



p-l 

E 

A=0 



K+M 



d m e (aAm!) 



Once again, by the Holder inequality 

2k-l 



'p-l 
vA=0 

Using (|TT|) . we obtain 

'p-i 



m=K+l 



p-l \ 2k - 2 p-l 

A=0 / A=0 



2k 



2k-l 



,A=0 
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We now have 

K+M 



£ 

A=0 



$ m e (aAm!) 



m=iC+l 

p-l K+M k 2k 

= ^2 E Htfm„e(aAraJ) J~| (-oAm,!) 

A=0mi,...,m2i c =Jf+li'=l u=k+l 

K+M k 2k p-l / / k 2k 

= e n aA i>* ! - e 

m 1 ,...,m 2 k=K+lu=l v=k+\ A=0 \ \f=l v=k+\ 

< P J k (K, M). 

Using Theorem Q we obtain 

\W a (K,M;L,N)\ 2M < pM 2k{e ~ l) N 2l{k ~ l) J e (L, N)J k (K, M) 

« p(MN) 2M N- f 'I^M- k ^ k+l \ 

and the desired bound follows. □ 

For example, for every fixed e > 0, choosing sufficiently large k and £, 
Theorem 0] yields a nontrivial bound whenever NM > p 1+£ . 

For K = L = 0, N = M = p — 1, choosing k = I = 2, we obtain that the 
bound of Theorem 0] is of the form 0(p 2 ~ l l 2A ). 

It is immediate that Theorem 0] combined with the Erdos-Turan relation 
between the discrepancy and the appropriate exponential sums (see |H IT5]) 
gives essentially the same the bound (with only an extra factor logp) on the 
discrepancy of the sequence of fractional parts 



m\n\ 



P 



K + 1 < m < K + M, L + l<n<L + N. 



As we have remarked, this improves in several directions a similar result 
from [7j. 

We also remark that Theorem |3] can be reformulated as an upper bound 
on the average value of sums W a (K, M; L, N) over a = 0, . . . ,p — 1. 
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4 Asymptotic Formulas for the Number of 
Solutions of Mixed Congruences with Fac- 
torials 

Let T r (K, M; L, N; A) denote the number of solutions of the congruence 

r 

7ij!mj! = A (mod p), 

i=i 

K + 1 < mi, . . . , m r < K + M, L + 1 < n h . . . , n r < L + N. 

Theorem 5. Let K, L, M and N be integers with 0<K<K + M<p 
and 0<L<L + N<p. For any positive integers k, t, r and s such that 
s < r/2 and N 2 > M > N 1 / 2 , we have 

(MN) r 
p 

•C ][fr-l+l/2s-(r-2s)/2e(k+l) jyr-l/2(s+V)-{r-2s)/2k{l+l) (r-2s)/2kl 



T r (K, M; L,N;X) 



Proof. Using the standard principle, we express T r (K, M; L, N; A) via expo- 
nential sums. Then 



1 

T r (K, M; L, iV; A) = - ^ (W a (K, M; L, N)) r e(-aA). 



p-i 

P a=0 

Separating the term (MN) r jp corresponding to a = 0, we obtain 



T r (K,M;L,N;X) 



(MN) r 1 J 



P 



<-Y^\W a {K,M-L,N)\' 



Pa=l 



1 P_1 

< max \W a (K,M;L,N)\ r - 2s -J2\W a (K,M;L,N)\ 2s 

Ka<»-1 r> 



i<p-l p 

a=0 

< F S (K,M;L,N) max \W a (K, M; L, N)\ 

l<a<p— 1 



r-2s 



Using Theorem El and Theorem EJ we finish the proof. □ 

Taking M = N, r = 7, s = 2, k = £ = 2 in Theorem we derive that 
any residue class A modulo p has N^p" 1 (+0(iV~ 17 / 12 p n / 8 )) representations 
in a form 

mibi! + rr^lr^! + . . . + m-j\nj\ = A (mod p), 
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with K + 1 < mi, rii, . . . , m 7 , n 7 < K + N (provided 0<K<K + N<p). 
In particular, each A is represented in the above form for K = and some 
N of the size N = (p 33 / 34 ) . 

For integers r > 0, A, L and iV with 0<L<L + iV<pwe denote by 
Q r (K, M; L, N; A) the number of solutions of the congruence 



mini + rij! = A (mod p), 



:i2) 



Ki<r 



with K + 1 < m < K + M and L + 1 < n, m, . . . , < £ + N. 

Theorem 6. Let K, L, M and N be integers with 0<K<K+M<p and 
0<L<L + N<p. Then for any integers k, I > 1, the inequality 



Q r (K, M; L, N; A) 



MN 



r+1 



P 



^ j^l-l/2e{k+l)j^r+l/2(ri+l)+l/2(r 2 +l)-l/2k(£+l)pl/2kl 

holds, where r± = \r/2\, r 2 = \_{r + 1)/2J. 

Proof. We express Q r (K, M; L, N; A) via exponential sums: 

p-i 

P 



1 

Q r (K, M; L,N;X) = -J2 W a (K, M; L, N) (S a (L, N)) r e(-a\). 



a=0 



Selecting the main term MN r+l /p, corresponding to a = 0, we obtain 

MN r+1 

Q r (K,M;L,N;X) 

V 



1 

« max | W a (K, M; L, N) | - V | S a (L, N) \ r . 

l<a<p— 1 p 

a=0 



Using Theorem H] and the same arguments as in the proof of Corollary |21 we 
conclude the proof. □ 

Taking M = N,r = 49, k = £ = 2 in Theorem |UJ we derive that any 
residue class A modulo p has iV 51 p -1 (l + O(A^ _4397/,3900 p 9//8 )) representations 
in the form 

mini + rii! + . . . + n i9 l = A (mod p), 
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with K + 1 < m, n, n x , . . . , < K + N (provided 0<K<K + N<p). In 
particular, each A is represented in the above form for K = and some N of 
the size iV = 0(p 8775 / 8794 ). 

One can also derive that for any e > there exists an integer r, such 
that any residue class A modulo p can be representated in the form (|12j) with 
1 < m < p e and 1 < n, n\, . . . , n r < p. 

We now combine TheorerrQ with the estimate (P) from our work and 
apply the method Karatsuba [5] of solving multiplicative ternary problems. 

For integers k, £, r > 0, A, L and iV with 0<L<L + N<p, we denote 
by Rk,£, r (K, L, S; M, N, T; A) the number of solutions of the congruence 

(mi! + . . . + rrik\)(ni\ + . . . + n^tj . . . t r \ = A (mod p), 

with K + 1 < m 1 ,...,m k < K + M, L + 1 < n 1 ,...,n e < L + N and 
S+ 1 <h,...,t r < S + T. 

Theorem 7. Let K, L, M, N , S, and T, A be integers with < K < 
K + M <p,0<L< L + N < p andO < S < S + T < p and A ^ (modp). 
Then for any integers k,£,r > 1, the following bound holds: 



Proof. Let X be the set of multiplicative characters modulo p, see jTSj. We 
have an analogue of (HJ) 



R kAr {K,L,S;M,N,T;X) 



^ M k-l/2+l/2(k+l) N £-l/2+l/2(£+l) T 3r/4 p r/8^ ogp y/4_ 




(mod p), 
(mod p). 



(13) 
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Therefore, we have 



R kAr (K,L,S;M, N,T;X) 

K+M L+N S+T 

= E E E 

mi,...,m(j=if+l rti,...,n^=L+l ti ,...,t r =5+l 



P 



^ X (K! + • • • + m k \)(m\ + ... + n e \)ti\ . . . t r \X x ) 
P xe* \t=S+l / 

if+M L+7V 

2j x( m i ! + • • • + m k \) ^ x( n i- + ■ ■ ■ + n e ] -)- 

mi,...,mf 1 =K+l ni,...,ni=L+l 

Separating the term M k N e T r / (p— 1) corresponding to the principal character 
Xo and then using the bound for the sum over t, we obtain 



R kAr (K,L,S; M,N,T;X) 



M k N l T r 



p — 1 

<T 3r/y/8 (logp) r/4_L_ y 
P — 1 ^— ' 



X6* 
Xt^XO 



K+M 

2j x(mi! + . . . + m fc !) 

mi,...,m fc =if+l 



L+JV 



^ x(ni! + ... + n^!) 

ni ,...,n^=L+l 



We see, from (|T3j) . that 

^tE 



xe* 



xCmil + .-.+mfe! 

mi,...,mi l ,=K+l 



< J k (K,M) 



(the above estimate is almost an equality were it not for neglecting the terms 
with m\\ + . . . + m k \ = (mod p)). Similarly, 



V 



xe* 



L+N 



^2 xi^l + . . . + n e \) 

n\ ,...,ri£=L+l 



< Je{L,N). 
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Therefore, by the Cauchy inequality, we see that 



Rk,£,r(K. 




« T 3p /y/ 8 (bgp) r / 4 (J k (K, M)J t (L, N)) 



1/2 



Using Theorem Q we finish the proof. 



□ 



Taking M = N = T,k = £ = 2 and r = 3 in Theorem we have that 
any A ^ (mod p) has N 7 (p- l)" 1 (l + 0(iV~ 17/1 V 1/8 (logp) 3/4 )) represen- 
tations of the form 



with K+l<m 1 , m 2 , m, n 2 , ti, t 2 ,t 3 < K + N (provided 0<K<K + N< 
p). The above asymptotic formula is nontrivial for any fixed e > and 

N > p 33/34+ £ _ 

One can also take k — 1, £ — 2, r = 4in Theorem [?l and obtain an 
asymptotic formula for the number of representations of the form 



with K + 1 < ni, ri2, ti, t 2 , ^3, tji, ts < K + N, which becomes nontrivial for 
N > p 18 / 19 + £ . However, our next result provides a stronger bound. 

Let us define R^ r (L, S; N, T; A) as the number of solutions of the congru- 
ence 



with L + 1 < rii, . . . , rig < L + N and S + 1 < ti, . . . , t r < S + T. That is, 
R ttr {L, S; N, T; A) = R QAr (0, L, S; 1, N, T; A). 

Theorem 8. Let K , L, M , N , S , and T , A be integers with < K < 
K + M <p,0<L< L + N < p andO < S < S + T < p and A ^ (mod p). 
Then for any integers £, r > 1 and < s < r, the following bound holds: 



(mj + m 2 !)(ni! + n 2 \)ti\t 2 \t^\ = A (mod p) 



(m! + n 2 ^t 1 \t 2 \t 2l \U\t 5 \ = A (mod p), 



(14) 



(n 1 \ + ... + n t \)t 1 \...t r \ = X (modp) 



R e , r (L,S;N,T;X) 



N e T r 



p — 1 



^-1/2+1/2(^+1) y(3r+s)/4-l/2+2 



—s-1 



> (r '- s)/8 (logp) (r - s)/4 . 



P 
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Proof. As in the proof of Theorem |HJ we derive 
R e , r (L,S;N,T;X) 



V 



1 / S+T \ ' L+N 

^T,^- 1 ) Exw E 

X €X \t=S+l / ni ,...,n e =L+l 



X{ni\ + . . . + n e \ 



Separating the term N e T r /(p — 1) corresponding to the principal character 
Xo, and also then using the bound (0), we obtain 



R i>r {L,S; N,T;X) 



p — 1 



«E 

X+XO 



S+T 



E x(t) 



t=S+l 



L+N 

E x(nil + ■ ■ ■ + n £ \ 

ni,...,ri£=L+l 



< T 3(r - s)/4 p (r_s)/8 (logp) (r " s)/4 



x¥=xo 



S+T 



E X(t) 



(=5+1 



E + • • • + nt\) 

ni,...,ri£=L+l 

Applying the Cauchy inequality, as in the proof of Theorem |Hl we derive 

_/V^T r 



R £ir (L,S;N,T;X) 



p — 1 



< T 3(r - s)/4 p (r " s)/8 (log p) (r " s)/4 {I S {S, T) J t (L, N) ) 1/2 , 

and using (j3J) together with Theorem Q we finish the proof. □ 

Taking N = T,r = 5,£ = 2 and s = 2 in Theorem |H1 we obtain that any 
A ^ (mod p) has N 7 (p - l)" 1 (l + 0{N- 35 / 2A p u ^(logp) 3 ^)) representa- 
tions of the form with K + 1 < rii,n 2 , ti, £2, £3, £4, £5 < K + N (provided 
0<iC<ir + iV"<p), which becomes nontrivial for iV > p 33 / 35+e . 

5 Remarks 



As we have remarked, a more careful examination of the function n(e) would 
lead to a substantial relaxation of the condition N 2 > M > N 1 ^ 2 of Theo- 
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rem El This however, does not affect the most interesting "diagonal" case 
M = N. 

Our method can also be used, without any substantial changes, to study 
the distribution of the products 

n 

Y[ f(j) (mod p), L + l <n<L + N, 

3=1 

where f(j) is rational function (welldefined modulo p for j = 1, . . . ,p — 1). 
In particular, with f(j) = j^ 1 one can estimate exponential sums and the 
number of solutions of some congruences with ratios of factorials n\/m\. 

Probably the most challenging open questions are obtaining a nontriv- 
ial upper bound on the exponential sums S a (L,N), and also obtaining an 
asymptotic formula (or at least a lower bound) on the number of solutions 
of the Waring-type congruence with factorials 

ri\\ + . . . + np\ = A (mod p), 

where L + l < ni, . . . , rig < L + N. Even in the case ofL = 0, N = p — 1 
these questions are still unsolved. Theorem |U1 seems to be the closest known 
"approximation" to a full solution of the Waring problem with factorials 
modulo p. 

As we have mentioned, our ability to obtain any extensive computational 
evidences is very limited. So, we dare not make any conjectures about pos- 
sible answers to the above questions. 
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